Kondo effect in triple quantum dots 
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Numerical analysis of the simplest odd-numbered system of coupled quantum dots reveals an 
interplay between magnetic ordering, charge fluctuations and the tendency of itinerant electrons 
in the leads to screen magnetic moments. The transition from local-moment to molecular-orbital 
behavior is visible in the evolution of correlation functions as the inter-dot coupling is increased. 
Resulting novel Kondo phases are presented in a phase diagram which can be sampled by measuring 
the zero-bias conductance. We discuss the origin of the even-odd effects by comparing with the 
double quantum dot. 
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PACS numbers: 72.15.Qm, 73.63.Kv, 72.10.Fk 

One of the main notions in molecular physics is the 
formation of extended molecular orbitals when atoms 
are brought together to bind into a molecule. Us- 
ing atomic manipulation capabilities of modern low- 
temperature scanning tunneling microscopes, the delo- 
calization of electrons can now be directly studied by 
adding individual atoms to a molecule-like chain of atoms 
on a fiat substrate Alternatively, one can study the 
transition from localized to extended behavior in the sys- 
tems of coupled quantum dots, which can also be consid- 
ered as a variety of artificial molecules. Unlike in atomic 
chains, bonding in coupled dots can be continuously ad- 
justed using the pinch gate electrodes |2j. 

The Kondo effect is a many-body phenomenon due to 
interaction between a localized spin and free electrons. 
It leads to an increased conductance through nanostruc- 
tures and it was first observed in single quantum dots 
(QD) 0- The conductance as a function of temperature, 
gate and bias voltages is in agreement with theoretical 
predictions that such dots behave rather universally as 
single magnetic impurities [J] and can be modelled us- 
ing single impurity Anderson and Kondo models. By 
analogy, the Kondo effect in multiple-dot systems can be 
described using multiple-impurity Anderson model. 

Measurements of the transport properties of multiple 
QD systems indicate that the conductance is significantly 
different for even and odd number of dots |5| . Double QD 
(DQD) systems are the simplest representation of the 
two-impurity Anderson model. Their phase diagram fea- 
tures single- impurity Kondo effect, two-electron Kondo 
effect and a suppressed-conductance state due to anti- 
ferromagnetic (AFM) coupling between the dots 0, • 
Good understanding of these different regimes was ob- 
tained using various numerical technia ues |sL 1^. Hfl Hi) . 

Much less is known about odd-number systems of cou- 
pled QD. Even the properties of the triple quantum dot 
(TQD), the simplest member of the family, have not 
been systematically resolved to date. There is clearly 
a regime of Coulomb blockade at intermediate tempera- 



tures 0, El • Some properties of the Kondo regime at 
low temperatures are known from previous studies using 
the perturbation theory [14| . the numerical renormaliza- 
tion group 0, scaling techniques Hr3l . the slave-boson 
mean-field theory 
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and embedding |l8j . Some of these 
methods lead to conflicting results for the same type of 
TQD, so further studies using reliable numerical tech- 
niques are required. 

In this Letter, we study the TQD using two com- 
plementary numerical methods: the constrained path 
Monte Carlo method (CPMC) and the Gunnarsson- 
Schonhammer variational method (GS) |2Cj. The con- 
ductance is calculated using the sine method |2l|. We 
analyze the low-temperature phases of the TQD embed- 
ded in series between two metallic leads, Fig. ^ 
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Figure 1: (Color online) The triple quantum dot embedded 
between two leads 

The TQD is modelled as a cluster of three equal QD 
described by Hamiltonians Hi = J2a e d n i<? + Un^riii, 
i = 1,2,3, where the on-site energy t c i can be varied 
using the gate voltage and U is the electron-electron (e- 
e) repulsion. The dots are interconnected with hopping 
matrix element t" and are symmetrically coupled to the 
leads with hopping matrix element t'. We choose the 
chemical potential in the middle of the band, so that the 
model is particle-hole (p-h) symmetric for + U/2 = 0. 

For large t" /U, the system is in the molecular- orbital 
Kondo regime when occupancy (n) — ^2 ia (ni a ) is odd. 
In the U — > limit, the conductance through the TQD 
as a function of the gate voltage exhibits three peaks cor- 
responding to resonant tunneling through distinct nonin- 
teracting levels: non-bonding, bonding and anti-bonding 
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molecular- orbitals with energies £±,0 — ±v2i") and hy- 
bridizations r ±)0 = |r, |r, where T = 2t' 2 /t. When 
interactions are switched on, the middle peak remains 
at + U/2 = 0, while the side-peaks are symmetrically 
shifted to e d + U/2 ~ ±{y/2t" + ^U) + <D{U 2 /t") in 
the molecular limit r — * 0. The [/-term is a conse- 
quence of the inter-orbital repulsion. The intra-orbital 
e-e repulsion is then given by U± ~ |i7 + 0(U 2 /t") and 
U ~ \U+ 0[U 3 /{t") 2 ]. The Kondo temperatures T ±j0 
for levels (±, 0) can be determined from the Haldane for- 
mula |22j for the single QD using the corresponding pa- 
rameters U± t o- 

The zero-temperature phase diagram of TQD exhibits 
several phases with enhanced conductance, Fig. [21 lu 
the molecular-orbital Kondo regime (shaded regions la- 
belled 'Ml' and 'M3') the conductance is G/G Q > 0.5 for 
low temperatures T <C T± j0j while for intermediate tem- 
peratures X±.o <C T -C T the system is in the Coulomb 
blockade regime with G ~ except along the border 
lines. Lightly shaded stripes of width ~ T/2 represent 
the transition regions, where G/Gq < 0.5. Phase 'M2' is 
a non-conductive even-occupancy spin-zero phase, where 
two electrons occupy the same molecular-orbital. 
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Figure 2: (Color online) Ml, M3: molecular-orbital Kondo 
regime with (n) ~ 1, 3. M2: non- conductive even-occupancy 
state. L3: local Kondo regime with (n) ~ 3. TSK: two-stage 
Kondo regime. K2: local double Kondo regime with (n) ~ 2. 

In Fig. |3 we present the zero-bias conductance G along 
with the total occupancy (n) as a function of the gate 
voltage for a range of t" and for fixed U/t = 1 and Y/t = 
0.18, calculated with various methods as discussed below. 
For t"/U > 0.2, the system is in the molecular-orbital 
regime. As the occupancy monotonically decreases from 
6 (full TQD) to (empty TQD), the conductance exhibits 
well resolved peaks when occupancy is odd and valleys 
when occupancy is even 0]. 

As the TQD coupled to the leads is Fermi liq- 
uid ^J, the zero temperature linear conductance is 
given by the sine formula derived in Ref. 0, G — 
G sin 2 [^(E -E n )]. Here G = 2e 2 /h and are 



the ground state energies of a large AT-site auxiliary ring 
with embedded interacting system, with periodic and 
anti-periodic boundary conditions, respectively. Ground 
state properties are determined using CPMC and GS 
methods. 

In the CPMC method [1J| , the ground state wave func- 
tion I^q) is projected from a known trial function \^t) 
using a branching random walk that generates an over- 
complete space of Slater determinants \(j>) and can be 
written as |\Po) = c ^\4>)i where > 0. To completely 
specify |*o), only determinants satisfying (^T\<t>) > are 
needed, because l^o) resides in either of two degenerate 
halves of the Slater determinant space separated by a 
nodal plane. In this manner, the minus-sign problem is 
alleviated. Extensive testing has demonstrated a signif- 
icant insensitivity of the results to reasonable choices of 
|*t> EEi. The CPMC calculations were performed 
on a ring of 100-180 sites. As the number of sites with 
interaction is small, CPMC produces ground state en- 
ergies with excellent precision, typically of the order of 
AE/E = 1Q- 6 . 
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Figure 3: (Color online) G/Go and (n) vs. ed + U/2 (propor- 
tional to the gate voltage) for various inter-dot hopping t" . 
Note enhanced energy scale in panels (d) and (e). Error bars 
of CPMC are smaller than the size of circles. For comparison, 
we also show the Hartree-Fock (HF) as well as U = results. 

In the GS method, the ground state energies are de- 
termined as the minimum energy of the variational wave 
function |W) = 5Za=i c ^-^ |6), where Pa are the pro- 
jection operators for the three sites of the molecule, e.g., 
p l = IIi,<r (!-"*,<*)) p 2 = n 1A (l-n hl ) IL=2,3^( 1 ~ n *< CT )' 
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. . . , Pq± = Yli a n i.<y- The vacuum state |5) corresponds 
to the Fermi sea of a non-interacting Hamiltonian H with 
renormalized parameters. The variational wave function 
is exact in the U — > and t' — > limits and the test re- 
sults for a single Anderson site perfectly match the exact 
Bethe Ansatz results [2ll 123. |25| . 

For t"/U > 0.2 the conductance obtained with CPMC 
and GS methods shows good agreement. For lower t" /U 
the CPMC method is no more applicable since due to the 
computational restriction on the system size its energy 
resolution is insufficient to describe the small Kondo scale 

m. 

The local regime emerges when t" /U decreases below 
~ 0.2. The orbital description breaks down as the or- 
bitals start to overlap which leads to qualitative changes 
in the system properties. The relevant interaction be- 
comes U rather than U± : o. As presented in Fig. Efd,e), 
G/Go ~ 1 near the symmetric point €d + U/2 = 0, while 
in the charge transfer region, |e<j + U/2\ ~ U /2, the con- 
ductance exhibits humps separated by dips. Due to larger 
U, the Kondo scale in the local regime decreases. Note 
that for t"/U < 0.2 the HF calculation starts to deviate, 
signaling the onset of the strong coupling regime. For 
even lower t" /U the HF curves differ qualitatively from 
correct results. 

The relation between the DQD and the TQD in the 
local regime is first considered for (n) ~ 2. This corre- 
sponds to the p-h symmetric point of the DQD, where 
the two relevant energy scales are known 6] to be the 
exchange energy J = 4(t") 2 /U and the Kondo tempera- 
ture Tik = T K (U,Ti), with Ti = T/2, which is the hy- 
bridization of a single impurity coupled to one lead [22j . 
In TQD, two electrons can reduce their energy by hop- 
ping within the TQD, therefore instead of J the relevant 
scale is the kinetic energy linear in t". The transition 
between the phase 'M2' to the phase 'K2' (see Fig. J2J) 
occurs when y/2t" ~ 2T 1K . The phase 'K2' is a double 
Kondo phase, where the spin of each electron is screened 
by electrons from the nearest lead. 

The behavior of the TQD near the p-h symmetric 
point, (n) ~ 3, is different from the behavior of the DQD 
near its p-h symmetric point, (n) ~ 2, due to the differ- 
ent properties of integer and half-integer spin states. It 
is also different from the behavior of the DQD near its 
charge-transfer points, (n) ~ 1,3, which exhibit Kondo 
effect only for large t" /U. The TQD is fully conducting 
at the p-h symmetric point for any t" and as t" /U is 
reduced the system goes continuously through three dif- 
ferent Kondo regimes. The molecular-orbital 'M3' regime 
has already been discussed. For low t"/U there are two 
different regimes of local Kondo physics, (i) In the lo- 
cal Kondo 'L3' regime for J > 2T\k 1 the TQD is an- 
tiferromagnctically rigid and electrons form an ordered 
chain with spin 1/2 that exhibits the usual Kondo ef- 
fect. Transition between molecular-orbital 'M3' and lo- 



cal 'L3' is quite soft and determined by the competition 
between kinetic it") and magnetic (J) scale, (ii) In the 
two-stage Kondo 'TSK' regime for J < 2Tik, each side 
QD (i = 1,3) couples with Ti to the nearest lead and 
their spin is screened at T < T\k, while due to very 
weak coupling the remaining spin is screened at much 
lower temperature T2K oc Tik exp(— cT\k/ J) with c ~ 1 

na- 
in Table [Q we give quantitative relations for lines be- 
tween different phases of the system. Boundaries between 
phases Ml, M2, M3, and L3 are determined in the molec- 
ular F — > limit to the lowest non-trivial order in U (or 
t"). The expressions are in excellent agreement with nu- 
merical results. We note that all transitions are smooth 
(cross-overs) and there are no abrupt phase transitions. 

Phase 1 Phase 2 Condition 



empty 


Ml 


Ml 


M2 


M2 


M3, L3 


M3 


L3 


M2 


K2 


L3 


TSK 



< v 
t" > u 



S~U/2 + s/2t" 

'5~ U/2 + 3{t") 2 /U, 

6~U/8 + V2t", 

'<J~ U/2-^/2t" + 3{t") 2 /U, t"<U 

J~U{l/4+^(U/t") 2 ), t">U 
V2t" ~ J 

V2t" + 3(t") 2 /U ~ 2min(Tix,ri) 
J ~ 2T 1K 



Table I: Definitions of boundaries and cross-over regions be- 
tween various phases in Fig. Here 5 = ea + U/2 and 
J = A(t") 2 /U. 

The finger-print of the Kondo physics are not only the 
broadened conductance peaks, but also the formation of a 
local moment in the TQD, which is then screened by the 
electrons in the leads. In Fig. 0] we present (S 2 ), where 
S z = J2 i=1 Si Z is the z-component of the total spin of 
the TQD (left panels), and local i = j, nearest neigh- 
bor {if) = (12) and next-nearest neighbor {if) = (13) 
spin-spin correlations (Si Z Sj Z ) (right panels). For large 
t" , formation of local moment (spin-doublet) is evident, 
but the spin is not saturated, 4(Sf) < 1. For t" — > and 
(n) ~ 2, increased 4(5^) > 1 illustrates formation of the 
'K2' phase with large local moment and the absence of 
inter-site spin correlations. The two electrons are rela- 
tively independent and therefore total spin is enhanced 
(note that for two uncorrelated spins 4:(S Z ) — 2). Each 
electron is coupled by the Kondo mechanism to the ad- 
jacent lead. For (n) ~ 3, we see that neighboring spins 
tend to anti-align while left-most and right-most spins 
align as t"/U is reduced, which represents formation of 
an antiferromagnetically ordered chain when we enter the 
'L3' phase. 

Charge fluctuations An 2 — ({n — (n)) 2 ) and the cor- 
responding charge susceptibility \c = — (nT/4)d(n) /dea 
exhibit double peaks at positions of the conductance 
peaks as an additional characteristic of the Kondo ef- 
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Figure 4: (Color online) Moments (Sf) vs. e d + (7/2, left 
panels (al)-(el). Spin-spin correlation functions (Si Z Sj Z ) for 
(ij) = (11), (22), (12), and (13) as calculated from the GS 
wave function, right panels (a2)-(e2). 



feet |2jj and start to build a dip in the symmetric point 
for t"/U < 0.2 (not shown here). This additionally sig- 
nals the transition from orbital to local Kondo regime 
and is in agreement with local spin formation shown in 
Fig.|H;d,e). 

In summary, we have determined the phase diagram of 
the TQD. For strong inter-dot coupling, the system be- 
haves as an artificial molecule. The extended "molecular 
orbitals" are filled consecutively and the system exhibits 
the usual Kondo effect when the number of confined elec- 
trons is odd. For weak inter-dot coupling, local spin be- 
havior is observed. The cross-over from extended to lo- 
cal Kondo physics is illustrated by the smooth evolution 
of the spin-spin correlation functions. In the vicinity of 
the p-h symmetric point, the TQD tends to an antiferro- 
magnetically ordered state when the inter-dot coupling 
is decreased. For extremely weak coupling, the Kondo 
screening of local moments occurs in two stages. In the 
charge transfer regime the left and right sites tend to form 
two relatively independent Kondo correlated states. 

We expect that in the limit of extremely weak coupling, 
other complex phases might exist in the immediate vicin- 
ity of the ed — point, where occupancy of the central 
dot changes abruptly. Furthermore, wc believe that the 
properties of all n dot systems (with n odd) are similar 
at the p-h symmetric point: as the inter-dot coupling is 
decreased, the molecular-orbital Kondo regime Mn is fol- 
lowed by an antiferromagnetically locked Kondo regime 
Ln and then by a two-stage Kondo regime. 
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